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Abstract 



An attempt to determine theoretically the highly non-linear current-voltage {I-V) 
characteristics of polycrystalline semiconductors, such as ZnO-based varistors, is made 
from the electrical properties of individual grain boundaries under dc bias. The role 
played by the fluctuations of double Schottky barrier heights at grain interfaces on 
driving electrical breakdown phenomena of macroscopic samples is pointed out in terms 
of a binary mixture model. An alternative trial form for the double Schottky barrier 
height is introduced to reproduce the breakdown voltage as well as the high non-linear 
coefficient a, where I ocV a . 



PACS numbers: 77.20, 77.50, 85.30 



A grain boundary in polycrystalline semiconductors (PCS) becomes electrically ac- 
tive as a result of charge trapping by gap states localized in the interface between adjacent 
grains. Such interface states may be created by the impurities, dislocations or interfacial 
defects and give rise to the appearance of double Schottky barriers M. The electrical con- 
duction through a grain boundary Schottky barrier may be determined by different con- 
ducting mechanisms, such as tunneling, field emission, space-charge limit, thermoionic 
emission, etc., (see, e.g., M), depending on the material and preparation. In (com- 
mercial) ZnO-b&sed varistors at room temperature, thermoionic emission is believed to 
be the dominant mechanism for electrical conduction through a (grain-boundary-grain) 
GBG junction [|, iH. 

Nevertheless, a theoretical analysis of the non-linear current-voltage (I-V) charac- 
teristics of a macroscopic sample of PCS from the point of view of electrical properties 
of individual GBG junctions is not trivial. This non-linearity is due to the complex 
charge and current feedback mechanisms occuring at grain boundaries which lead to a 
large variety of potential barriers. Thus, the barrier height (and, in correspondence, 
the breakdown voltage) of a GBG junction differs from one to another as shown by 
microelectrode measurements ||, [Tj. The random nature in size and position of grains 
is revealed by electron micrograph experiments |2|, |3|, |6], [7|, |8| . 

The simplest model proposed to characterize ZnO-varistors is due to Levinson and 
Philipp [[J. In their model all ZnO-GBG junctions are assumed to be identical and 
regularly arranged in a cubic lattice. To a crude approximation, the I-V characteristics of 
PCS then became completely determinated by the properties of one GBG junction. More 
recent studies, that included parallel or series GBG structures or used the percolation 
approximation, have also neglected fluctuations of double Schottky barrier heights |S|, [X0f| . 
Thus, the question of how the non-linear coefficient a, where / oc V a , and the electrical 
breakdown voltage of PCS may be affected by the distribution of barrier heights at grain 
boundaries is still open. 

In this work we made the first attempt to determine in simple terms the highly non- 
linear I-V characteristics of ZnO-b&sed varistors under dc bias by considering the non- 
linear electric characteristics of single grain boundaries. For that purpose, we use a 3D 
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effective medium approximation (EM A) [[II], [L2| for a binary mixture of grain boundary 
Schottky barriers, having different heights, in conjuction with the thermoionic emission 
mechanism (TEM) for the carrier transport through grain boundaries. An alternative 
trial form for the voltage dependence of a single barrier height is also introduced. Herein, 
we study as a function of barrier probability the variations in the macroscopic (i) non- 
linear coefficient a and (ii) electric breakdown voltage Ve of a PCS sample featured by 
the binary mixture. 

The quantity a is commonly defined as 

a(V) = ^r . (1) 
K ' d(lnV) K ) 

For the breakdown voltage Vb, however, there are somewhat different definitions. For 
the present purposes, it is natural to define Vb as the voltage related to the maximum 
value of ct(V), i.e. as the solution of 

/ da . . 

( W )W.= ■ (2) 

Accordingly, by a = o:(Vb) we shall denote the maximum value of the non-linear coeffi- 
cient. 

In the traditional TEM the current that is injected through a potential barrier is 

B § § B 1 

I = A*T 2 e [ fe s T J (l-e , (3) 

where A* is Richardson's constant, V is voltage drop across a junction, Ef is the Fermi 
energy in the grain bulk (counted from the bottom of the conducting band), T is the 
temperature, e is the elementary charge and 4> is the potential barrier height that depends 
on V. 

Solving the Poisson equation, within the simplest framework of the double Schottky 
barrier model (i.e. neglecting effects due to holes and interface processes), such that the 
interface charge is independent of V and time, <fi(V) is found to be || |13| 

pV 

4>{V) = O (1 - — f ; (eV<40 o ) , (4) 
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where <p = <p(V = 0) is the potential barrier height at zero applied V. The voltage 
V c = that is related to 4>(V C ) = 0, is often interpreted as the critical voltage for 
electric breakdown in a single GBG junction ||. 

Nevertheless, according to definitions ([I]) and (0) and to the I-V characteristics of 
an individual GBG junction, i.e. Eqs.(|3]) and (f|), the breakdown voltage to a first 
approximation can be estimated from 

V B ^ , (5) 

and, in correspondence, the non-linear coefficient then becomes 

oc ~ — ; — ^ • (6) 
2k B T v ; 

Thus, the value of Vb given in Eq.fl5|) is two times less than the above V c value. The 
measured breakdown voltages are really about V c |5], || |9|. For a typical value of <f) ~ 
0.75 eV at room temperature, Eq.(^[) leads to a ~ 15, which is not high enough to 
explain the observed strong non-linearity of I-V characteristics of typical ZnO-b&sed 
varistors (a ~ 25 to 50 or more) 0. 

In Fig.l (left-hand side axis) the behaviour of <j)(V) as given by Eq.(^j) is compared 
with measured data normalized to a single junction (0 O = 0.73 eV; T = 300 K) JT3[ . 
This figure shows the poor agreement obtained with experimental results which have 
only a slight voltage dependence at small values of applied V and suddenly decrease 
in the vicinity of the breakdown region. In this figure the principal features of Eq.(^) 
using Eq.([|) are also depicted for the same set of parameters (right-hand side axis). 
We note that at -j^j, < < 1, i.e. the prebreakdown region, the I-V characteristics 
obtained using (|3]) and @ present an ohmic behaviour with activation energy equal 
to (j> . Without paying much attention to this ohmic region, the curves in Fig.l have 
been drawn to illustrate electrical breakdown phenomena at a GBG junction. In fact, at 
y ~ t hi s fi 

gure shows (solid line, right-hand side axis) that the current in the TEM 
suddenly increases as V increases in accord to Eq.@. The value of a obtained from 
the curve /(V) is exactly -^^f, (c.f. Eq.fl^)). For V > it is experimentally found 
that the current reaches the value limited by the resistance of the ZnO grains. The 
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I-V characteristics are then described by a second ohmic law which is not accounted by 
Eq.(§. 

The above discussion concerns an individual GBG junction only. As mentioned there 
exists many junctions in PCS which are not (necessarily) identical and have different 
o 's. So, we determine next by numerical approach the macroscopic I-V characteristics 
of PCS. To achieve this we consider a 3D-EMA for a binary mixture of GBG junctions 
with non-linear I-V characteristics and use the TEM for d.c. conduction through single 
GBG junctions. 

The EMA, originally proposed in Ref . [Ill , is a self-consistent scheme used to deter- 
mine the conductivity of inhomogeneous materials (see also [0). ^ consists on replacing 
the random net by a homogeneous lattice, where all lattice bonds have the same con- 
ductance J2- Choosing the lattice spacing to be unity, then J2 is also the macroscopic 
conductivity of the homogeneous material which can be derived self-consistently as dis- 
cussed in Ref. fTT|, |12[ . However, it is important to mention here that this approach has 
been initially applied to linear systems with constant bond conductances a\, 02- Noting 
that in this work a\ and 02 of GBG junctions depend on the applied voltage V, it is 
straighforward to repeat the EMA procedure to get a system of self-consistent equations 
for determining the sample conductivity as a function of V. In this case, the predic- 
tion of the EMA for J2 of a 3D binary distribution of bond conductances is obtained 
from 

P » 1 (P) + 2E 1 W <r 2 (l") + 2£ ' ( '> 

and 

vr = 3E+^)-^(v) {i = l2) (8) 

where p is the probability of having bonds with conductance o\. 

We represent next the ZnO-grains as lattice sites and the ZnO-GBQ junctions as 
bonds. Applying Eqs.(|3|) and @ to obtain o\ and 02, each of which is derived assuming 
different values for and (j)^ , we can then estimate the conductivity of the system as 
a function of probability p by solving self-consistenly Eqs.(^) and (§). 
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Figure 2 shows our results for the binary mixture conductivity. These are obtained 
from EMA within the probability range 0.1 < p < 0.9 by using typical data for ZnO- 

,(1) ,(2) 

based varistors, i.e. = 9 and |^ = 30 J/J. In the event of an extremely narrow 
distribution of one type of barrier heights (either p=0 or p=l), the predictions of the 
EMA can be shown to be asymtotically correct [0. To this end we emphasize again 
that we are dealing with a binary mixture of grain junctions having positive p < 1. 

From Fig. 2 we get interesting results for the sample conductivity over a wide p- 
interval. For p < 0.3, the binary mixture conductivity presents a behaviour similar to 
the case of single GBG junctions with highest double Schottky barriers (p^ (c.f. Fig.l). 
On the other hand, when p > 0.5, the sample is essentially determined by the lowest 
double Schottky barriers having (f>£\ In the intermediate region 0.3 < p < 0.7, as 
discussed below, it is worthwhile to mention that there is an assymetric contribution to 
the sample electrical properties from each GBG component. 

As it is well known, at the EMA percolation threshold value for the component one, 
i.e. p^ ~ 1/3 ||12|| , the binary sample conductivity slightly deviates from a smooth 
increasing on varying V. Similarly, at the percolation threshold value for the second 
component, i.e. p^ ~ 1 — |, the sample also deviates from a smooth behaviour on 
increasing V. This is due to deficiences in the EMA which is known to give unaccurate 
predictions at these particular junction probabilities |I2J. Furthermore it must be noted 
that, for all p- values considered, the sample conductivity continues to be an increas- 
ing function of applied voltage until -as we shall see next- electric breakdown occurs. 
For voltages V > Vb, J2 tends to the grain conductivity which is chosen as unit for 
conductivities in our calculations. 

Having seen from Fig. 2 that, to some extent, the EMA is able to reproduce the typical 
non-linear I-V characteristics of ZnO-based varistors, it is then possible to give an 
estimate of Vb and a for the macroscopic samples. In Figs. 3(A) and 3(B) we summarize 
the p-dependences of Vb and a for the binary mixture, respectively. The full lines 

,(2)_ ,(1) 

correspond to different (quality factor) x = ( i) of barrier heights; namely (1) x — 2, 

,(i) 

(2) x — 1 an d (3) x = 1/2, where Y~f = ^ is fixed throughout calculations. The three 
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curves displayed in each figure are respectively normalized to the electrical breakdown 
voltage -j^f = i^f = ^0 anc ^ to ^ ne non-linear coefficient = -^h^ = 5 due to 
component one. Since we have chosen cp^ > from Eqs.(|5|) and (|6|) we then have 
> VfP and > At p = 0, the sample Vb and a become identical to those 
values derived from a junction with barrier height At p ^ 0, it is clear that the 
sample Vb and a varies between vjp — > Vg and 

As visualized from Figs. 3 (A) and 3(B), the sample Vs{p) and a(p) exhibit sharp 
variations around p = 1/2, especially so at x > 1. When % << 1, a smooth dependence 
on p of Vb and a is found on decreasing the probability of a W -component. In the present 

(2) 

cases the sample a reaches a maximum value of a/ 2 -* ~ 2k B T ' ^ * s i m P or t an t to note that, 
when the conductances (or, alternatively, (f>)P and (j)^) of the components of the binary 
mixture have nearly same probability p ~ 0.5, there is an asymmetric contribution 
of the components to the non-linear coefficient and electric breakdown voltage of the 
sample. This asymmetry becomes large on increasing \- These results imply that the 
components play different roles on driving macroscopic I-V characteristics. 

The component with barrier height (f>^ becomes important only at p < 0.5 and 
dominant for p < p£\ On the other hand, for all p > 0.5, the sample a and Vb 
are enterely determined by the component with lower barrier height, i.e. <p^\ U P to 
an applied voltage V such that the second barrier conductivity can not be neglected. 
Simultaneous contributions from both components at higher V determine the behaviour 
of sample (/-V) characteristics in the range p > 0.7. Such assymetric features of the 
electrical properties of PCS follow from the monotonic behaviour of the I-V curves shown 
in Fig. 2. Moreover, the sample a in Fig. 3 is dominated by the potential barriers with 
average or lower heights. The higher double Schottky barriers do not play considerable 
role in this case. 

In the experiments by Olsson and Dunlop ]7), several types of junctions between 
ZnO grains were distinguished. The various types of junctions characterized by differ- 
ent intergrain structures revealed different I-V behaviour. The most frequent interface 
structures were found to be of "type A", i.e. ZnO boundaries containing a thin inter- 
granular amorphous .Bj-rich film, "type B", i.e. boundaries between ZnO grains not 
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involving any second intergranular phase but containing Bi atoms, and "C type", i.e. 
ZnO grains separated by an integranular region containing Bi 2 Oz. Since types A and 
B exhibited similar I-V characteristics, the present theoretical GBG junctions with (f>jp 
may be related to the "C-type" junctions whereas the GBG junctions with (f)^ to the 
"B-type" junctions (by neglecting possible differences in the bias dependence for these 
types of GBG). Further motivation to approximate macroscopic PCS as a binary mix- 
ture of grain boundary Schottky barriers is obtained from the experiments by Tao et 
al. ||. Our theoretical results in this case may represent the large majority of grain-to- 
grain junctions characterized from an electrical point of view as "good" junctions and 
"bad" junctions. As can be seen in Fig. 3, our work predicts that for large proportions of 
"good" GBG, i.e. no less than 70%, the sample non-linear I-V characteristics become 
completely determined by such junctions. 

The presence of free electrons in the conduction band, the hole current (created by hot 
electrons at high enough applied voltages), and interface processes (electron captures, 
emission and hole-electron recombination at interface states) (see e.g. f|, [L4], [0|), 
might be important to derive a more satisfactory expression for the barrier height 4>(V) 
than Eq.([|). The new expression should lead to obtain higher values of Vb and a for 
GBG junctions, and from it, for the ZnO-based varistors too. This complicated problem 
has been solved numerically for a single GBG junction [^, |T5| . Up to now, there is not an 
analytical expression for 4>(V) which, on the one hand, gives desired results for a single 
junction problem [IT]] and, on the other hand, could be basic for studing macroscospic 
non- linear I-V characteristics of ZnO-based varistors as it concerns us here. 

As a simple example, let us consider the alternative trial form 

0(V)=0 o e /32 ^(l-/? 1 ^) 2 ; (PieV<4> ) . (9) 

It can be seen in Fig.l by dotted lines (left-hand side axis) that this choice (though not 
unique) allows us to fit better experimental data [13| than using Eq.(|4]) from a simple 
adjustment of two parameters, namely f3\ and Clearly, the first term of an expansion 
series of the exponential function in Eq.(|9]) leads to Eq.(^) for low V . The following 
terms in such a series might represent those additional effects that become important at 
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the high voltage region p 2 eV ~ (f> . Theoretical I-V characteristics derived using Eq.([|) 
for the thermoionic emission current and Eq.(^]) for 4>(V) are also illustrated in Fig.l by 
dotted lines (right-hand side axis) for completness. 

If we take Pi and P% such that, at << 1, the barrier heigh becomes weakly 
dependent on V, (i.e. J^- | v ^ ~ 0). Then, we immediately obtain p 2 = 2/3i and Eq.(Q) 
for 4>(V) becomes dependent only on Pi which may be chosen to correlate the voltage 
V B (~ with a high non-linear coefficient a for the GBG junctions. 

In fact, from definitions (|l|) and (0), for enough high barriers we have 



Vi 



B 



V2Pi 



e 



a * 32p 2 1 (l-4Pi)e^^ f . (10) 

For Pi fa 0.19 we get the desired Vb and a fa 1.28^^ which is more than twice higher 
than a from Eq.(|6|). The coefficient a becomes independent of Pi, but depends on 
the zero- voltage barrier heigh <f> . If, as before, <p D fa 0.75 eV then Eq.([H]) gives a ~ 
35 for a single GBG junction at room temperature. Hence, we feel that Eq.(|9|) may 
be a convenient expression to mimic 4>(V) in these systems in the absence of simpler 
approaches to include effects due to holes and interface processes. 

Let us focus again briefly on the EMA. In Figs. 3 (A) and 3(B) (dotted lines) we 
plot new values of the sample Vb{p) and a(p) by using Eqs.@ together with Eq.(f|) and 
X = 0.5 for a comparison with the solid lines 3. We immediately note that the behaviour 
of the sample a and Vb is the same than in the previous case by applying Eq.(f|). But, 
even more important, according to relations in Eq.(|T0"D a and Vb reach higher values. 
The same conclusion can be drawn for different values of x ( n °t shown). 

To conclude some comments are in order. We have analysed the non-trivial problem 
of the non-linear I-V characteristics of a macroscopic sample of PCS from the point of 
view of electrical properties of individual GBG junctions. In real percolation systems 
a network of non-linear elements at breakdown is known to be filamentary. In order 
to take into account of this effect, a more sophisticated model which includes random 
fluctuations of (multiple) barrier heights must be considered. In this work we have 
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restricted ourselves to a simple binary mixture model of PCS (ZnO)-grain interfaces 
and used the TEM for the carrier transport through the grain boundaries. The use of 
a binary mixture approach being motivated by the experimental findings of Refs.|7|. |8|]. 
In the case of a narrow distribution in fluctuations of barrier heights ^ [2(J], the I-V 
characteristics of varistors may became very similar to the single junction data and the 
simplest model of Levinson and Philipp suffices || . 

We have found that features in the ^/-dependence of 4> and fluctuations of dou- 
ble Schottky barrier heights at (rather than the electrical conduction through) grain 
boundaries play an essential role to achieving highly, non-linear I-V characteristics of 
macroscopic PCS samples. For this reason, we have introduced a simple trial form for 
the voltage dependence of a single barrier height in the absence of simpler approaches 
to mimic effects due to holes and interface processes. We have assumed that the bias 
dependence of GBG is the same for both components by considering either Eq.([|) or 
Eq.(^|). As shown in Fig.3, different 4>{V) functions will clearly lead to obtain different 
results depending on the chosen form. The 3D-EMA for non-linear I-V characteristics 
of PCS as introduced in this work presents an asymmetric contribution of barrier heights 
to the sample a and Vb around p = 0.5. This asymmetry becoming large on increasing 
the difference in the barrier height values of the components. 

These interesting findings may be useful guidelines to simulation studies of electric 
breakdown to a network of random potential barriers with non-linear I-V characteristics 
|IR |TT| I n addition, similar EMA calculations may be also extended to 2D systems 
such as polycrystalline semiconductor films |J. In the latter, we shall expect analogous 
results because 2D-EMA equations are equivalent to those in Eqs.(|7|) and (|8|) for 3D 
samples but with a small change of coefficients. 
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Figure 1: Double-Schottky barrier height at, and dc current through, a single 

grain boundary as a function of applied voltage V for <p = 0.73 eV at T = 300 K. 

_ JLl 

Here I a = A*T 2 e k s T . Full lines are obtained from Eqs.(|^) and (|]) and dotted 
lines from Eqs.(|3]) and (|9[). Squares represent experimental data for <f>(V) p~3|j . 

Figure 2: Results for the binary mixture conductivity as a function of applied 
voltage V for several probabilities p of having bonds with conductance o\. Typical 

,(1) ,(2) 

data for ZnO-based varistors are used: = ® anc ^ ftr = 

Figure 3: (A) Normalized voltage Vb at which electrical breakdown occurs and 
(B) normalized non-linear coefficient a for a binary mixture as a function of prob- 
ability p of having component one. In these plots (1) \ = 2, (2) \ = 1 and (3) 
X = 1/2. Full lines are obtained using Eqs.(0) and (|J). Dotted lines are obtained 
using Eqs.(|3|) and @ with x — 1/2. 
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